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Abstract 

In a recent experiment at the vuv free-electron laser facility at DESY in Hamburg, the generation 
of multiply charged ions in a gas of atomic xenon was observed. This paper develops a theoretical 
description of multiphoton ionization of xenon and its ions. The numerical results lend support to 
the view that the experimental observation may be interpreted in terms of the nonlinear absorp- 
tion of several vuv photons. The method rests on the Hartree-Fock-Slater independent-particle 
model. The multiphoton physics is treated within a Floquet scheme. The continuum problem 
of the photoelectron is solved using a complex absorbing potential. Rate equations for the ionic 
populations are integrated to take into account the temporal structure of the individual vuv laser 
pulses. The effect of the spatial profile of the free-electron laser beam on the distribution of xenon 
charge states is included. An Auger-type many-electron mechanism may play a role in the vuv 
multiphoton ionization of xenon ions. 

PACS numbers: 32.80.Rm, 31.15.-p, 87.50. Gi 
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I. INTRODUCTION 



An electron bound to an atom experiences electric forces, which on average point toward 
the atomic nucleus. If the atom is placed in a static electric field, the electronic states 
become unstable, as the potential arising from the superposition of the atomic and the 
external electric field enables electron emission via tunneling. This picture basically remains 



lation period 




valid even if the external electric field is oscillating, at least as long as the osci 
of the electric field is long in comparison to the electron tunneling timescale 
For a typical valence-electron binding energy of the order of 10 eV and a laser photon energy 
of about 1 eV — i.e., in the ir regime — the tunneling picture is meaningful for intensities of 
10^'^ W/cm^ or higher. However, at high radiation frequency (or at low intensity), this 
picture fails, and it is more appropriate to adopt a multiphoton picture 0, 3, S An 
atomic electron can be ejected following the absorption of a discrete number of photons. 

The development of free-electron laser (FEL) llOll facilities at several places in the world 
Q has spurred substantial theoretical interest Q 

until very recently, no radiation sources have been available to experimentally probe strong- 
field physics at vuv or shorter wavele ngth s. The situation changed when the VUV-FEL 
at DESY, Hamburg, began operation [3; l^^- In one of the first experiments with this 
exceptional radiation source, xenon clusters were exposed to the intense vuv laser pulses 
and were observed to absorb a surprisingly large number of v uv p hotons per atom 2^, a 
finding which was explained in terms of inverse bremsstrahlung . No evidence for atomic 
multiphoton processes was found in these first measurements. Measurements using a more 
sensitive detector, however, revealed the creation of multiply charged ions even in a gas of 



free, unclustered xenon atoms [26[ . An experimental time-of- flight mass spectrum, averaged 
over 100 consecutive VUV-FEL pulses, is shown in Fig. ^ S^]. The photon energy in that 
experiment was 12.7 eV with a peak intensity of approximately 10^^ W/cm^. 

In this paper, we present a theoretical description of the interaction of xenon atoms 
with an intense pulse of vuv photons. Our findings lend support to the interpretation 



26[, is a consequence of 



that the occurrence of multiply charged ions, as reported in Ref. 
vuv multiphoton ionization — a scientific first at a photon energy of more than 10 eV. We 
describe our computational method in Sec. |H] In Sec. IIIIl the results of our calculations 



are presented and compared with experiment. We conclude with Sec. IIVI Atomic units are 
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used throughout, unless otherwise noted. 



II. COMPUTATIONAL METHOD 



Using an independent-electron model, Geltman 



(see also Refs. 



29|,|30|) was able to 



arrive at a satisfactory representation of the ionic charge distributions observed in several 



intense- laser experiments at photon energies of 6.4 eV and below 
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33, 



34, sj- 



Within 



the framework of independent particles, each electron moves in the field of the atomic nucleus 
and in a mean-field generated by the other electrons. The best such mean-field derives from 

n 

the Hartree-Fock self-consistent-field method !36]. However, the Hartree-Fock mean-field is 
nonlocal, due to the exchange interaction, and therefore cumbersome to work with. Slater 
js?! introduced a local approximation to electron exchange, which is the principle underlying 
the well-known Xa method 38^ . 

In this work, we use the Hartree-Fock-Slater code written by Herman and Skillman j^^, 
which in the past has proved advantageous for atomic photoionization studies (see, for 
example, Ref. |40|). The resulting one-electron potential, Vhs('^)) is a central potential 
(even for open-shell systems), which satisfies 



Vhs('^) -Z/r , r ^ 



(1) 
(2) 



for an atom of nuclear charge +Z and overall charge +q. In the original program of Herman 
and Skillman, the Xa parameter was set to unity, in accordance with Ref. j37|. We adjusted 



that parameter slightly (see Table I!]), so that the ca^ 



culated bp binding energy in Xe, Xe"*", 



Xe"'""'", and Xe^"*" agrees with experimental data 4l|, |42, |43|, |4^. The situation is a little 
more complex for Xe^"^ and Xe^"*". The experimental ionization potentials of these ions, 
as determined from measured electron- impact ionization thresholds j3, li^ vary by 



as much as 18 eV. Utilizing the complete-active-space se 
implemented in the ah initio program package MOLPRO 



f-consistent-field (CASSCF) code 

nn 

48| , I49l | ; an active space consisting 



of the 4(i, 5s, 5p, 6s, and 6p orbitals; and an effective core potential 



we obtained 



ionization potentials of 53 and 64 eV for Xe^"*" and Xe^+, respectively. The value for Xe^+ 

n n 

agrees with Ref. |45| ; the one for Xe^"*" is in agreement with Ref. j47| . As indicated in Table 
EJ we base the Xa parameters for Xe^"*" and Xe^+ on the ionization potentials quoted above. 
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It should be mentioned that while our CASSCF calculation includes scalar relativistic effects 
in an approximate manner (through the effective core potential), spin-orbit splittings have 
not been taken into account. 

We represent the radial one-electron Hamiltonian 



H 



(I) 

AT 



1 d^ + , , 



(3) 



in a finite-element basis set jsi, 52, S^ls^Q^H' which is described in detail in Ref. j57[. In 
the calculations reported here, 4801 finite-element basis functions were employed, spanning 
a radial grid from rmin = to rmax = 80 bohr radii. For each orbital angular momentum 
quantum number, /, considered, the first 50 eigenfunctions Un,i{r) of H^J^}^ were computed 
that have eigenenergies en,i at least as high as the energy £5^0 of the 5s level. More strongly 
bound levels are assumed to be fully occupied by core electrons and are not considered any 
further. The calculated eigenfunctions M„,i(r) satisfy the boundary conditions Mn,i(^min) = 
and 'Un,i(rmax) = 0. At the current level of approximation, atomic multiplet structure is 
absent. All electrons that can be associated with the quantum numbers n and I have the 
same energy, 6n,i, within our model. In particular, there is no energy dependence on the 
magnetic quantum number, m. 



In order to treat the problem of electron emission, a complex absorbing potential (CAP) 

irjW^r), is added to 
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60, ml ^ 
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701, 
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complex absorb] 

HHQ, -i 



the atomic Hamiltonian. The real, positive parameter rj is the CAP strength. The local 
one-electron potential W{r) is chosen here as 



W{r) 




(4) 



where c = 4 a.u. in this paper (which places the absorbing potential right outside the ionic 
core). The CAP absorbs the outgoing electron and renders the associated wave function 
square-integrable. Given a complete basis, there exists for a resonance eigenstate (a Gamow 
vector) of the physical Hamiltonian with Siegert eigenvalue i^res = -E'r ~ ir/2, an eigenvalue 
E{r]) of the CAP-augmented Hamiltonian such that lim^^o+ E(r]) = Erpf, l62l. In a finite 
basis set, r^opt must be found, satisfying 

dE{r]) 

»?opt 



V- 



drj 



mmimum . 



(5) 
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E{riopt) provides an approximation to the Siegert energy -Eres, from which the resonance posi- 
tion, E-£i, and the resonance width, F, can be extracted. An improved strategy, which is used 
lere, consists in analyzing the function E{ri) — ridE{r])/dr] and minimizing \ri'^d'^E{'r])/dr)'^\ 



Within the framework of quantum electrodynamics the Hamiltonian describing the 
effective one-electron atom interacting with the electromagnetic field, in the presence of the 
CAP, reads 

H = Hat + Hem + H, - i7]W , (6) 



where 



i^AT = -^V^ + ^HSW , (7) 

-f^EM = ^ uja\yak,\ , (8) 

fc,A 



Hi = X -J^^yy^ {^k,xak,x - el ^^al^^^ . (9) 

fc,A 

Here, -f^EM represents the free electromagnetic field. Hi the interaction term in electric-dipole 
approximation (in length gauge). The operator al.^ {ci'k,\) creates (annihilates) a photon 
with wave vector k, polarization A, and energy uj = k/ a {a is the fine-structure constant). 
We use the symbol x for the atomic dipole operator. V in Eq. (jH)) denotes the normalization 
volume of the electromagnetic field, and Cfc a indicates the polarization vector of mode A;, A. 
Let be the number of photons in the laser mode, so that the laser intensity is given by 

intensity being measured in units of Jq = -Eh/(^oc^o) ~ 6.43641 x 10^^ W/cm^ [Ei^: Hartree 
energy; to: atomic unit of time; Qq: Bohr radius). We can now combine the atomic 
eigenstates, 'ipn,i,m = {un,i{r)/r)Yi^rn{^,'^), with the Fock states of the laser mode, \N — fi) 
(/i = 0, ±1, ±2, . . .), to form basis vectors \^n,i,m,fi) = \4'n,i,m) \N — jj). Assuming linear po- 
larization, the matrix representation of the Hamiltonian H [Eq. ©] in the basis {|$n,«,m,ft)} 
is diagonal with respect to m. It also has a rather sparse structure with respect to n, Z, and 
/i. The only nonzero matrix elements are 
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I Hat + Hem I'&n.,;,™,^*) 

{^n,l,m,^\ Hi |$n',Z',r?i,/i+l) 
(^n,Z,m,^+l I Hi I $n',/',m,,^) 
\^ n,l,m,fj.\ \^n' ,l,in,fj./ 



V2TraI {lpn,l,m\ Z \lpn\l',m) 
V27ral {lpn,l,m\ Z \i'n',l',m) 
{'4^n,l,in \W \tlJn',l,m) ■ 



(11) 

(12) 
(13) 
(14) 



The energy Nu of the unperturbed laser field has been subtracted from the right-hand side 
of Eq. (jllj) : the relatively high intensity {N ^ |/i|) has been exploited in the coupling matrix 
blocks (Eqs. (jl2j) and (|13p): and a unitary transformation has been applied that renders the 
matrix -ff at + -f^EM + iJi real s ym metric. The complete matrix H is complex symmetric 



and of the Floquet-type 



U real s yra meiric in 

H 22, H 29, so.lsi 



8al8l| (see, for example, Refs. |82l |8 



85 



and references therein for other computational approaches to atomic strong- field physics). 
In our calculations, fi runs from = to yUmax, the minimum number of photons needed 
to photoionize (referred to as N.P. in Table H}. Thus, since electron emission can only take 
place after the absorption of /imax photons, it is sufficient to apply the CAP only to the 
/imaxth diagonal block, i.e., ($n,«,m,/x| W \^n',i,m,fj) is set to for 7^ /i^ax- 



III. CALCULATIONS 



As a test of our method, we determined the one-photon ionization cross section of neutral 
Xe at a photon energy of 12.7 eV. In this calculation, s, p, and d waves were included, and 
/^max was set to 1. Let us first consider ionization of the 5p, m = level. After assembling 
the corresponding real symmetric matrix H^t + -f^EM + Hi, those 120 eigenvectors of this 
matrix were computed that have the largest overlap with the initial-state vector | $5,1,0,0)- 
The ?7-dependent complex symmetric eigenvalue problem of H [Eq. ©] was then solved in 
the subspace of the previously calculated eigenvectors of the real part of H. The chosen 
subspace size of 120 provided converged results and, at the same time, allowed for extremely 
fast optimization of the parameter 77 [see Eq. (0) and the text following it]. 

The rj trajectory of the resonance energy in the complex plane, for an intensity of 1 x 10^^ 
W/cm^, is shown in Fig. |21 On the basis of this graph, a dynamic Stark shift of 7.31 x 10^^ 
a.u. and an ionization rate of 75^,0 = 3.05 x 10^^ a.u. are found. Proceeding in a similar 
fashion, the ionization rate of the 5p, m = ±1 levels is calculated as 75p,±i = 2.01 x 10~^ 

6 



a.u. at 1 X 10^^ W/cm^, so that the m-averaged ionization rate 

75p = (75p -1 + 75p,o + 75p,i) /3 (15) 
is 2.35 X 10~^ a.u. We then calculate the total ionization rate as 

Tsp = 75p (6 - g) , (16) 

where, for neutral Xe, g = (g is the atomic charge). This procedure is approximately valid 
also for g > 0, since spin-orbit interaction ensures that the 6 — g 5p electrons are uniformly 
distributed over m = —1, 0, +1. 

We have checked that T^p is a linear function of the intensity, /, in the vicinity of 1 x 10^^ 
W/cm^. Thus, within the Herman- Skillman-based independent-particle model, the one- 
photon ionization cross section of neutral Xe, at a photon energy ^ of 12.7 eV, is 119 Mb. 
This result, which differs from the experimental cross section [jST,] by a little more than 
a factor of 2, has been confirmed by us using the same independent-particle model, but 
treating the continuum problem with an R-matrix code Q| (see also Ref. j3])- 

In order to test whether our CAP-Floquet program is also capable of describing multi- 
photon physics, we investigated the two-photon ionization cross section of neutral Xe at 6.42 
eV. Experimentally, this is known to be about 4 x 10"^'^ cm"^ s js^. With atomic s through 
/ waves and /imax = 2, we calculated a cross section of 6.2 x 10~^° cm^ s. Other calculations 
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of this quantity, which are similarly accurate, are reported in Refs. 

We calculated the (g + l)-photon ionization cross section, a^+i, of Xe''^ at a photon 
energy of 12.7 eV (see Table H]) following a strategy analogous to the one described above 
for neutral xenon: 

^2 (Xe+) = 4.5 X 10-^^cmS, (17) 

ag (Xe++) = 4.7 x 10"^^ cm^ s^ , (18) 

a4(Xe^+) = 1.8 X 10-^^^cmSs^ (19) 

as (Xe^+) = 1.1 X 10-^^^ cm^" s^ , (20) 

(Xe5+) = 2.4 x 10"^^^ cm^^ ^5 _ (21) 

The spatial profile of the VUV-FEL beam in Hamburg, perpendicular to the beam axis, 
has a Gaussian shape [2^. Let {p,ip,z) denote cylindrical coordinates with respect to that 
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axis. The intensity near the focus (at z = 0) may then be written as 

41n2 / 41n2 ^\ , , 

where 

A(^) = A^l + (z/zoy (23) 

is the 2;-dependent full-width-at-half-maximum of the Gaussian beam profile. In the exper- 
iment described in Ref. A = 20 fim. The beam divergence was 17 mrad [2^, from 
which we estimate that zq = 1.2 mm. The time-dependent pulse power is represented in Eq. 

m by Pit). 

Since we are interested in the nonlinear response of Xe ions to the vuv laser pulses, it is 
not permissible to use for P{t) the pulse shape obtained after averaging over many pulses. 
The temporal shape of the individual FEL pulses has not been measured so far, but reliable 
simulations of the FEL performance exist which are able to reproduce measured 

^EL parameters and which, in addition, provide information about temporal pulse shapes 
9^ . Ten representative, simulated pulses are shown in Fig. |31 9^. We see that while the 
averaged pulse may appear approximately Gaussian (with a pulse width of about 50 fs), the 
individual pulses are not. 

InCud.. ..e„u..o„ facto, of 0.2 Q, wh.h fa., .to accou„f fKe B.fe .Becf.v., 
of the mirrors used to focus the FEL beam into the xenon gas, we solve, for each of the 
pulses shown in Fig. 13 the rate equations 



no{p,z,t) = -CTi no{p,z,t), 

ni{p,z,t) = ai no{p,z,t) - a2[ ni{p,z,t), (24) 



UJ \ iU ^ 

2 / r/ J.N \ 3 



n2{p,z,t) = 0-2 ( 1 ni{p,z,t) - as \ ) n2{p,z,t), 



for the probabilities nq{p,z,t) of finding Xe'^"'" at time t and position [p, z) [ip arbitrary]. 
(The thermal motion of the ions on a timescale of 100 fs may, of course, be neglected.) The 
initial conditions are no(p, z,t —00) = 1 and nq{p, z,t —00) = for g > 0. Let k, stand 
for the gas density in the interaction region. Then the total number of Xe''"'" generated by a 
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given laser pulse reads 

Nq = 2tck, / dz dp pnq{p, z,t ^ +00) . (25) 
In the laser experiment at DESY, k, = 2.8 x 10^^ atoms/cm^, Zmin = ~1 nim, and ^max = +1 



mm 22 1- 



We calculated Ng for each of the 10 laser pulses in Fig. |21and then determined the average 
number of Xe'^"'" generated per laser pulse, < Ng >, which is depicted in Fig. ^^It is difficult 
to assess whether this result can already explain the measurements in Ref. ^6|. According 
to the calculation, Xe"*" and Xe"*"*" dominate by far, which is consistent with the fact that 
the detector response to these two ions appeared to be saturated in the experiment 3|. 
The mass spectra in Ref. j2^ have not been calibrated to account for the specific detector 
response to ions of different charge states [2^, so they may not be linearly related to the 
actual ion production rate. Moreover, several of the experimental parameters we employed 
in our calculation are, in fact, not known very precisely. Among these are the gas density, 
K, and the spatial beam width, A 27]. It should be mentioned, in addition, that the FEL 
was not operating under optimal conditions when the data in Ref. j3| were taken |27|. 
Therefore, the laser pulses our calculation is based on (Fig. ^ are, on average, more intense 
than in the experiment. 

The theoretical model we use also suffers from shortcomings. As mentioned earlier, mul- 
tiplet splittings of the valence shell are not considered, which means the intermediate bound 
states influencing the multiphoton ionization cross sections may not be sufficiently accurate. 
The spectral width, Alu/uj, of the laser pulses of, on average, 1 % |27] — which is broader 
than the Fourier limit of a 50-fs pulse by an order of magnitude — is not included in the 
present treatment. In the rate equations, Eq. (j21I), excited-state populations, phase effects, 
as well as nonsequential multiphoton processes are neglected. The latter, however, may be 
expected to be strongly suppressed in the vuv regime, as confirmed by the measurements in 



Ref. 



26| 



Before concluding, we would like to mention an interesting many-electron effect that leads 
to an enhancement of the multiphoton ionization rate in the xenon ions. In Xe"*", it requires 
11.3 eV to excite one of the two 5s electrons to the 5p shell. (Within the Herman- 
Skillman model we find 10.9 eV.) One can, therefore, envisage, as one of the paths leading 
to two-photon ionization of Xe"*", the excitation of a 5s electron by a ffist vuv photon and 



the subsequent excitation of one of the 5p electrons to a virtual state bound in the channel 
associated with the 5s hole, as illustrated in Fig. Due to electron correlation, one of 
the remaining 5p electrons can fill the 5s hole, and the excited electron is ejected into the 
continuum. This is a kind of Auger decay of the inner-valence excited ion, resulting in the 
formation of Xe++. A similar scenario is conceivable for the more highly charged xenon ions. 

The contribution of this to the multiphoton ionization cross section of Xe'^"'" can be crudely 
estimated as follows. We choose | $5,0,0,0) as initial-state vector in our Floquet code, but 
instead of utilizing a CAP, we simply assign an autoionization width Fauto (i-e. we add 
— irauto/2) to those diagonal elements in the /imaxth diagonal block that satisfy + £5,1 — 
£5,0 > 0. This condition implies that the energy released in the 5p 5s transition is 
sufficient to transfer the electron with quantum numbers n and / into the continuum. A 
typical lifetime of an inner-shell hole is of the order of 10 fs (or shorter), so we set the 
autoionization width to 65 meV. 

Let us call 75^ the 5s-mediated ionization rate determined in this way. The total 5s- 
mediated ionization rate of Xe"^^ is then, approximately, 

T5s = 2^5skr-q) ■ (26) 
6 

The factor of 2 in this expression is needed since there are two 5s electrons in all Xe ions 
considered here. If we assume that the six 5p spin orbitals have equal probability of being 
occupied, then the probability that the 5p, m = spin orbital with the right spin is un- 
occupied is q/6. After the virtual excitation of one of the 5s electrons, there are 7 — q 5p 
electrons available for the absorption of the remaining /imax — 1 photons. 

We add the cross sections obtained in this way to the respective cross sections in Eqs. 
(I17|) - (j2ip . thereby neglecting interference effects. The results are 



^2 (Xe+) 


= 4.6 X 10" 


-49 4 

cm s , 


(27) 


as (Xe++) 


= 2.0 X 10" 




(28) 


(74 (Xe3+) 


= 3.3 X 10" 


"^^^ cm^ s^ , 


(29) 


a, (Xe^+) 


= 3.7 X 10" 


"i4^cmi0s\ 


(30) 


ae (Xe^+) 


= 6.4 X 10" 


"i^9cmi2s\ 


(31) 



The quantities most significantly affected by the 5s-mediated ionization mechanism are the 
three-photon ionization cross section of Xe+"^ and the five-photon ionization cross section of 
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The ionic charge distribution computed using the multiphoton ionization cross sections 
in Eqs. (P7 |) -(| ^ is displayed in Fig. IHl Now the production of Xe^+ is clearly visible even 
on a linear scale: About 3x10^ triply charged xenon ions are produced per laser pulse. The 
number of Xe^^ ions per laser pulse is more than 3000, which is also not particularly small. 



IV. CONCLUSION 



We have investigated in this paper multiphoton ionization of atomic xenon and its ions at 
a photon energy of 12.7 eV, a radiation intensity of order 10^^ W/cm^, and a pulse duration 
of about 50 fs. A recent experiment employing the VUV-FEL at DESY has demonstrated, 
under these laser conditions, the generation of xenon charge states of up to 6+ 2o|. In 
the infrared, even a pulse that is three orders of magnitude longer (with about the same 
intensity), produces charge states no higher than Xe^"*" 

Using an effective one-particle model, in combination with the Floquet concept and a 
complex absorbing potential, we have calculated vuv multiphoton ionization cross sections 
that refer to the absorption, by a 5p electron, of as many photons as are needed to ionize it. 
We have also estimated the influence of 5s excitation on the ionization cross sections and 
found that it may be significant. Although we grant that the model we applied is not ideal for 
describing many-electron phenomena (a better many-body calculation would be desirable), 
the result of our estimate indicates that at vuv photon energies multiphoton ionization is 
driven, at least partly, by electronic many-body physics. Focussing on the behavior of a 
single active electron does not appear to be sufficient. 



Taking a rate-equation approach and utilizing simulated FEL pulses j95| , we determined 
the average number of Xe''^ ions produced per vuv laser pulse. This step depends heavily 
on a number of important experimental parameters |22l]. In our calculation, thousands of 
Xe^"*" ions are found to be generated per pulse, and a correspondingly higher number for the 
lower charge states. Hence, the experimental observation of multiple ionization of xenon, 
Ref. appears compatible with the nonlinear absorption of several photons. We conclude 
that multiphoton physics is indeed relevant for some processes driven by the intense vuv 
beam of the free-electron laser in Hamburg. For more quantitative comparisons, it will be 
desirable for future experiments to obtain a calibrated ionic charge distribution, as well as 
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more detailed information about the FEL pulse properties. 
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FIG. 1: Experimental |27l | ion detection signal is plotted versus the ionic time-of- flight. Free- 
electron laser pulses with a peak intensity of about 10^^ W/cm^, a duration of approximately 100 
fs, and a photon energy of 12.7 eV produce, in an atomic xenon beam, ionic species of various 
charge states Q- This trme-of-flight mass spectrum ^ is an average over 100 consecutive FEL 
pulses. 
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FIG. 2: The dots in this figure correspond to a specific eigenvalue of the family of complex sym- 
metric matrices H{r]n) [Eq. ©], where r/„ = — 1)/(k — 1) [n = 0, ... ,99; 5 = 3 x 10~^; 
K = 1.1]. In the absence of photons and absorbing potential, this eigenvalue equals the 5p, m = 
level of atomic xenon (which defines the origin in the figure). Due to the interaction with 12.7-eV 
photons, this level is shifted as well as broadened. The point of stabilization of the r] trajectory 
implies a dynamic Stark shift of 7.31 x 10^^ a.u. and an ionization rate of 3.05 x 10^^ a.u., at a 
radiation intensity of 1 x 10^^ W/cm^. 
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FIG. 3: The power in the DESY free-electron laser pulse is shown as a function of time. The vuv 
laser pulses supplied by the free-electron laser source at DESY, Hamburg, are ultrashort, with an 
average width of about 50 fs, and intense, with a pulse energy of order 10 
shown in this figure [o^ have been calculated using an FEL simulation program 
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2c I . The pulses 
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FIG. 4: The average number, < A''^ >, of Xe'^^ ions produced per VUV-FEL laser pulse, as 
calculated on the basis of the multiphoton ionization cross sections in Eqs. H17 |) - (|21|) . the laser 
pulse properties [Eqs. ((^ . and Fig. the rate equations in Eq. (EU, and the integral over 
the interaction volume in Eq. ()25() . Note the logarithmic scale along the ordinate of the inset. 
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FIG. 5: Schematic depiction of energy levels associated with the excitation of xenon ions from 
the 5s subshell. For instance, Xe"*" can be photoionized via the absorption of two vuv photons by 
a 5p electron. A second ionization path, alluded to in this figure, involves a photon absorption 
that virtually excites a 5s electron to the bp shell. Simultaneously, a 5p electron is excited to an 
autoionizing state associated with a Rydberg series converging to the 5s threshold of Xe"'"''. 
Electron correlation then induces a transition from 5p to 5s accompanied by the emission of an 
electron. 
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FIG. 6: The same as in Fig. ^ the only difference being the fact that the multiphoton ionization 
cross sections underlying the calculation were taken from Eqs. ()27() - (|31|) . not ()17 |) -(|21 |) . < A'^g > is 
the average number of Xe'^"'" ions produced per VUV-FEL radiation pulse. 



TABLE I: Xa parameters employed to reproduce the 5p ionization potential (I. P.) of Xe^"*", i.e. 
the energy needed to remove one electron from Xe"^"^ and generate Xe^^"*"^^^ in its ground state. 
Also shown is the minimum number of 12.7-eV photons (N.P.) needed to ionize Xe"^"*". 
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